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1 Introduction 



The study of three-dimensional anharmonic oscillators has raised a considerable amount of interest 
due to its wide applications in molecular physics [1-4]. The Morse potential is commonly used for 
anharmonic oscillator. However, its wavefunction does not vanish at the origin. But those of Mie type 
and pseudo-harmonic potentials do. On the other hand, the Mie type potential possess the general 
features of the true interaction energy [1]. The Mie type and pseudo-harmonic are two special kinds of 
exactly solvable power-law and inverse power-law potentials other than the Coulombic and harmonic 
oscillator. The Hamiltonians with inverse power potential l/r n are studied in different fields. For 
instance, path integral, 1/N expansion and super-symmetrical approaches are also studied for this 
type of potentials [5-8]. Moreover, there are some applications of Mie type potentials such as equation 
of bulk metallic glass for Zr^NbgCuuNi^Beis in solid state physics [9], and other applications in 
solid state physics [10], molecular spectroscopy [11], fluid mechanics [12], the interatomic interaction 
potential in molecular physics [13], l/r n potential in three dimensions [14]. 

The Nikiforov-Uvarov (NU) method [15] is introduced for the solution of the hypergeometric 
type second order differential equations [16] appeared in the time-independent problems. Recently, 
NU-method has been used to solve Schrodinger equation for some well known potentials [17-22], 
Dirac and Klein-Gordon equations for the Coulomb and some exponential potentials. In the present 
work the Schrodinger equation is solved by the NU method for Mie potential with any value of 
angular momentum £. As an illustration, energy eigenvalues are computed for N 2 , CO, NO and CH 
molecules. 

The organization of the paper is as follows. In section II, NU method is briefly introduced. Mie 
potential calculations are done by using NU method in section III. Results are discussed in section IV. 

2 The Nikiforov-Uvarov Method 

The NU-method developed by Nikiforov and Uvarov (NU-method) is based on reducing the second 
order differential equations (ODEs) to a generalized equation of hypergeometric type [15]. It is 
provided us an analytic solution of Schrodinger equation for certain kind of potentials. It is based 



2 



on the solutions of general second order linear differential equation with orthogonal function [16]. 
The one-dimensional Schrodinger equation is reduced to a generalized equation of hypergeometric 
type with an appropriate coordinate transformation s = s(r). Then, the equation has the form, 



where A(s) = 44 and B(s) = J^r- In Eq. (1), a(s) and a(s) are polynomials at most second 
degree, and f(s) is a polynomial with at most first degree [15]. The wave function is constructed as a 
multiple of two independent parts, 



1>(s)=<f>(s)y(s), (2) 

and Eq. (1) becomes 



a{s)y"(s) + T(s)y'(s) + Xy(s) = 0, (3) 

where 



<j(s)=Tr{s)^-{ln<Ks)), (4) 



ds 



and 



A is defined as 



r(s) = f(s) + 2tt(s). (5) 



Xn + nr + H w - 1 ) g "] =0, n = 0,l,2,... (6) 



We determine 7r(s) and A by defining 



fc = A-7r'(s). (7) 
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From the Eqs. (4-5), n(s) becomes 



= (^) ± \] (^) 2 -& + k* (8) 

7r(s) has to be a polynomial of degree at most one in the Eq. (8), so the expression under the square 
root must be the square of a polynomial of first degree [15]. This is simply possible only if its 
discriminant is zero. After defining k, one can obtain n(s), r(s), <j>(s) and A. If we look at the Eq. (4) 
and the Rodrigues relation 



R <7 n 

where C n is normalizable constant and the weight function satisfy the relation as 



±[a(s)p(s)]=r(s)p(s). (10) 



where 



7T S 



<f>(s) a(s)' 



(11) 



3 Mie Potential Calculations 



The Mie-type potentials are given by 



V(r) = e 



k fa 



t (a\ k 



(12) 



— k \rj £ — k \r 

where e is the interaction energy between two atoms in a solid at x — a, and £ > k is always satisfied. 

We solve the one-dimensional Mie potential [1-2] with £ = 2k combination. By choosing the spe- 
cial case k — 1, corresponding to a Coulombic-type potential with an additional centrifugal potential 
barrier, we get the following form: 



V(r) = V'„ 



if")'- 

2 \r J \r 



V () = 2ek 



(13) 



where V is the dissociation energy and a is the positive constant which is strongly repulsive at shorter 
distances. The radial part of the Schrodinger equation for a diatomic molecule potential is 



!*LLJL ( 2 —\ £(?+l)h 2 
2/1 r 2 dr \ dr J 



2^r 2 



+ V(r) 



R n i(r) = E n£ R n£ (r), 



(14) 



where fi is the reduced mass of the diatomic molecules, n denotes the radial quantum number (n and 
£ are named as the vibration-rotation quantum numbers in molecular chemistry), r is the internuclear 
separation. Substituting the explicit form of V(r), we get 



d 2 R nE (r) | 2 dR nE (r) | 2/i 
dr 2 r dr % 2 



E„ 



nl 



v„ I ( a - 



2 Vr 



a 
r 



Rni{r) 



2/ir 



l)h 2 

J—R ne (r)=0. (15) 



By defining the following variables: 



E 



P = ~T2 V o a 
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F ( 2 ° a + 



2/i i 

'ft 5 

+ l)ft 2s 
2/x / 



The Schrodinger equation takes the simple form: 



(16) 
(17) 

(18) 



d 2 2d 1/on \ 

+ -— + -(eV-/3r- 7 

r rtr r z ^ > 



Rni{r) = 



(ir 2 r <ir r 2 

From Eq. (19), it is clear that a(r) = r, f = 2, a = e 2 r 2 — f3r — 7. We find 7r(r) as 



(19) 



1 

7r r) = — ± { 

v ; 2 



ier - |V 4 7 + 1, ^2 = 



-/3 + eV-i-47; 

-/3 - eV-1-47- 



and r(r) can be written as 



(20) 



r r 



1 + 2ier + V4 7 + 1, fci 
1 - 2ier + ^7 + 1, ^2 



-/3 + eV-1-47; 

-P - £y/-l -4 7 . 



(21) 



For appropriate solutions r'(r) < 0[15],weuse k 2 = —(3 — e^—l — Aq, n 2 = —\—%er — + 47 
and r(r) = 1 — 2ier + ^/l + 47 for obtaining eigenvalues and eigenf unctions. Using Eqs. (6-7) we 
find energy spectrum as 



*„ = -(^) : 2n+1 + 

Wavefunction is calculated from Eqs. (9-11). 



8/x /^o2 l(l + l)ft 2 
\| + h 2 { 2 + 2/i 



-2 



(22) 



0(r) = r 



i(-l+^/I+4^) e -i£r 



and 



(23) 



p(r) = r 

and the radial part of the wave function becomes 



y/l+A-y -2ier 



(24) 



= C n r^- 1 +^ I +^)e-^Lr^(2i£r) 
If we rewrite Eq. (22) in atomic units, we obtain 



(25) 



E n = -V 2 a 2 (2n + 1 + J (2£ + l) 2 + 2V a : 



(26) 



4 Conclusions 

We have investigated the analytical solution of the Schrodinger equation for Mie potential. By using 
a special case of Mie potential asm = 2k and k = 1, the problem is reduced to a Coulombic potential 
with the additional centrifugal barrier. Energy eigenvalues and corresponding eigenfunctions are 
obtained by the NU method, energy eigenvalues are computed for N 2 , CO, NO and CH molecules. 
They are listed in Table 1 . 
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Table 1: Calculated energy eigenvalues of the Mie potential for N 2 , CO, NO and CH diatomic 
molecules with different values of n and £ in eVThe data for potential parameters are taken from 



[22]. 



State (n) 


i 


N 2 


CO 


NO 


CH 








0.00271631 


0.0254063 


0. 


0205538 


0.041601; 


1 





0.0810235 


0.0756383 


0. 


0611503 


0.120556 




1 


0.0810268 


0.0758723 


0. 


0613568 


0.122185 


2 





0.13411 


0.125172 


0. 


101132 


0.194768 




1 


0.134351 


0.125403 


0. 


101335 


0.1963 




2 


0.134833 


0.125867 


0. 


101742 


0.199356 


3 





0.186486 


0.17402 


0. 


140511 


0.264609 




1 


0.186719 


0.174248 


0. 


140712 


0.266052 




2 


0.187194 


0.174704 


0. 


141112 


0.26893 




3 


0.187907 


0.175387 


0. 


141713 


0.273228 


4 





0.238152 


0.222196 


0. 


1793 


0.330417 




1 


0.238385 


0.222421 


0. 


179498 


0.331777 




2 


0.238854 


0.22287 


0. 


179892 


0.334491 




3 


0.239557 


0.223544 


0. 


180484 


0.338544 




4 


0.240495 


0.224442 


0. 


181272 


0.343916 


5 





0.28913 


0.269712 


0. 


217511 


0.392495 




1 


0.289361 


0.269933 


0. 


217705 


0.393779 




2 


0.289824 


0.270377 


0. 


218094 


0.396341 




3 


0.290518 


0.271041 


0. 


218677 


0.400166 




4 


0.291443 


0.271927 


0. 


219453 


0.405238 




5 


0.292599 


0.273036 


0. 


220424 


0.41153 
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